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I. INTRODUCTION 



Chiral gauge theories may display an anomalous breaking of Lorentz and CPT invariance 
if these theories are defined on a spacetime manifold M with nontrivial topology [1, 2]. One 
example of an appropriate manifold is M = M 3 x S 1 , as long as the chiral fermions have 
the correct boundary conditions (spin structure) over the compact space dimension. This 
so-called CPT anomaly has been established also for two-dimensional "chiral" U(l) gauge 
theories over the torus T 2 , where the Euclidean effective action is known exactly [3, 4]. 
The main points of the CPT anomaly have been reviewed in Ref. [5] and some background 
material can be found in Ref. [6]. 

The crucial ingredient of the CPT anomaly is the existence of a compact, separable space 
dimension with appropriate spin structure. For the case that this dimension is closed, there 
are corresponding noncontractible loops over the spacetime manifold, which must therefore 
be multiply connected. Up till now, attention has been focused on flat spacetime manifolds 
M with topology M 3 x 5 11 or I x T 3 and with Minkowski metric g^x) = rj^ and trivial 
vierbeins e£(x) = 5^ . Here, the noncontractible loops occur at the very largest scale. But 
there could also be noncontractible loops from nontrivial topology at the very smallest scale, 
possibly related to the so-called spacetime foam [7-14] . The main question is then whether 
or not a foam-like structure of spacetime could give rise to some kind of CPT anomaly. 

In the present article we give a positive answer to this question. That is, we establish the 
CPT anomaly for one of the simplest possible manifolds of this type, namely M = Rx M 3 , 
where the three-dimensional space manifold M 3 is flat Euclidean space M 3 with one straight 
line R removed. This particular manifold M 3 has a single "puncture" and arbitrarily small 
noncontractible loops. 

But it is, of course, an open question whether or not spacetime really has a foam-like 
structure and, if so, with which characteristics. Let us just consider one possibility. Suppose 
that the spacetime manifold M. is Lorentzian, Ai = M. x Ai^, and that the topology of 
is multiply connected and time-independent (here, time corresponds to the coordinate 
x° gR). The idea is that the multiply connectedness of .M3 is "hard- wired." The advantage 
of this restriction is that, for the moment, we do not have to deal with the contentious issue 
of topology change; cf. Refs. [12, 14]. 

Physically, we are then interested in the long-range effects (via the CPT anomaly) of 
this foam-like structure of space on the propagation of light. The exact calculation of these 
anomalous effects is, however, not feasible for a spacetime manifold M. with many punctures 
(or similar structures). We, therefore, introduce a model for the photon field over M 4 which 
incorporates the basic features found for the case of a single puncture. The model is relatively 
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simple and we can study the issue of photon propagation in the large- wavelength limit. 

The outline of this article is as follows. In Sec. II, we establish the CPT anomaly for a flat 
spacetime manifold, where the spatial hypersurfaces have a single puncture corresponding 
to a static linear "defect." (A similar result is obtained for a space manifold with a single 
wormhole, which corresponds to a static point-like defect. The details of this calculation are 
relegated to Appendix A 1.) 

In Sec. Ill, we present a model for the photon field which generalizes the result of a 
single puncture (or wormhole). We also give a corresponding model for a real scalar field. 
Both models involve a "random" background field over R 4 , denoted by go(x) for the scalar 
model and by g±(x) for the photon model. For the photon model in particular, the random 
background field gi(x) is believed to represent the effects of a static spacetime foam (a more 
or less realistic example for the case of permanent wormholes is given in Appendix A 2). In 
Sec. IV, we discuss some assumed properties of these random background fields. 

In Sec. V, we calculate the dispersion laws for the scalar and photon models of Sec. 
III. In Sec. VI, we use observations of gamma-rays from explosive extragalactic events to 
place an upper bound on the typical length scale of the random background field g\ (x) and, 
thereby, on the typical length scale of the postulated foam-like structure. In Sec. VII, we 
present some concluding remarks. 

II. CPT ANOMALY FROM A SINGLE DEFECT OF SPACE 
A. Chiral gauge theory and punctured manifold 

In this section, three-space is taken to be the punctured manifold M 3 = lx (R 2 \ {0}) = 
R 3 \ R. The considered three-space may be said to have a linear "defect," just as a type- 
II superconductor can have a single vortex line (magnetic flux tube). The corresponding 
four- dimensional spacetime manifold M = R x M 3 is orientable and has Minkowski metric 
= diag(l, — 1, — 1, — 1) and trivial vierbeins e a ^(x) = 8®. This particular spacetime 
manifold M is, of course, geodesically incomplete, but the affected geodesies constitute a 
set of measure zero. At the end of Sec. II B, we briefly discuss another manifold, N, which 
is both multiply connected and geodesically complete. 

The main interest here is in quantized fermion fields which propagate over the predeter- 
mined spacetime manifold M and which are coupled to a given classical gauge field. We will 
use cylindrical coordinates around the removed line in three-space and will rewrite the four- 
dimensional field theory as a three-dimensional one with infinitely many fermionic fields, at 
least for an appropriate choice of gauge field. This procedure is analogous to the one used 
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for the original CPT anomaly from the multiply connected manifold R 3 x S 1 . 

The gauge field is written as B^(x) = gB b (x)T b , with an implicit sum over b, where 
g is the gauge coupling constant and the T b are the anti-Hermitian generators of the Lie 
algebra, normalized by tr (T b T c ) — — \ 5 bc . For the matter fields, we take a single complex 
multiplet of left-handed Weyl fermions ip a (x). As a concrete case, we consider the 50 (10) 
gauge theory with left-handed Weyl fermions in the 16 representation. This particular chiral 
gauge theory includes, of course, the Standard Model with one family of quarks and leptons. 
Incidentally, the anomalous effects of this section do not occur for vectorlike theories such 
as Quantum Electrodynamics. 

In short, the theory considered has 

[G,R L ,M, el{x) ] = [ 50(10) , 16 , M 4 \ R 2 , ^ ] , (2.1) 

where G denotes the gauge group, Rl the representation of the left-handed Weyl fermions, 
M the spacetime manifold, and e a (x) the vierbeins at spacetime point x G M. The action 
for the fermionic fields reads 

/fcrmion = / d*X $ & K i (d, + B„) ^ , (2.2) 
J M 

with 

(d°* a ) = (l,-ff) , (2.3) 

in terms of the Pauli spin matrices a 1 , a 2 , a 3 . Natural units with h — c — 1 are used 
throughout, except when stated otherwise. 

Next, introduce cylindrical coordinates (p,(f>,z = x 3 ,t = x°), with 

x 1 = p cos , x 2 = p sin cj) , 

d f} = cos (p di + sin 82 , d<f, = —p sin d\ + p cos 82 , 

a p = cos a 1 + sin a 2 , = — sin cr 1 + cos <r 2 . (2.4) 

The two spin matrices of Eq. (2.4) and a z are explicitly 
/ — e~ i( ^ \ ( ie - ^' 






For the spinorial wave functions in terms of cylindrical coordinates, we take the antiperi- 
odic boundary condition 

i/>(p,<j> + 2Tr,z,t) = -i/>(p,<j>,z,t). (2.6) 
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An appropriate Ansatz for the fermion fields is then 

oo 

MpA,z,t)= {e +i(n - 1/2) *X< + W,*)^ +) +e +Kn+1/2) *X { n\p,z,t)i>i- ) }, (2.7a) 

n=— oo 

oo 

&(p,«M,f) = E { e_i(n_1/2) ^i +) (p^^)^ +) +e- i ^ +1 ^xi- ) (p,^t)^- ) }, (2.7b) 

n=— oo 

with constant spinors 

W +, ) = (j). W- ; ) = (°), ('^ +, ) = (o). (*'") = (!)• < 28 > 

and anticommuting fields Xn ± ' ) which depend only on the coordinates p, z, and t. The 
unrestricted fields x^\p-> z ->t) an d Xo^(P> z i t) of Eqs. (2.7a,b) will play an important role 
in the next subsection. 

B. CPT anomaly 

In order to demonstrate the existence of the CPT anomaly, it suffices to consider a special 
class of gauge fields (denoted by a prime) which are ^-independent and have vanishing 
components in the <fi direction; cf. Ref. [1]. Specifically, we consider in this subsection the 
following gauge fields: 

B' <p = B'^p,z,t) = 0, B' m = B' m (p,z,t), m = p,z,t. (2.9) 

Using the Ansatz (2.7) and integrating over the azimuthal angle 0, the action (2.2) can be 
written as 

;>oo poo poo 

/fermion =2vri £ / dp p / d Z / dt { ( d t + B[ ) X « + & ( % + 3 ) ^ 

- ( d p + b' p ) x h - ( 9 P + 5; ) x <+> - ( 9, + ) X « 

+ ( d z + B' Z ) X ^ + - p ((n~ 1/2) - (n + 1/2) } , (2.10) 

with e a positive infinitesimal. 

This action can be interpreted as a three-dimensional field theory with infinitely many 
Dirac fermions, labeled by n G Z . These three-dimensional Dirac fields are defined as 
follows: 

fj n = (v^xi +) , V^PXI'^I , (2.11) 



V^px^ ] \ 

V^PXn ] J 
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with 7-matrices 

= a\ V^ 3 , 7 2 = -i^\ (2.12) 

which obey 

7°7 1 7 2 = i, = 2rT, (2.13) 

for /i, v — 0, 1, 2, and [rf v ) = diag (1, —1, —1). The action (2.10) is then given by 

/fermion = i Y\ {fin i 7^ ( ^ + B' ) ?? n + ^ Vn Vn \ , (2. 14) 

Jm 3 I y J 

where t, p, z have been renamed y°, y 1 , y 2 , respectively, and B' has been defined as B' = B' t , 
B[ = B' B' 2 = B' z . The relevant three-dimensional spacetime manifold M 3 = Ex M. >0 x M 
has no boundary and is topologically equivalent to M 3 . 

The n = sector of the theory (2.14) describes a massless Dirac fermion in a background 
gauge field, whereas the n ^ sectors have additional position-dependent mass terms. At 
this moment, there is no need to specify the gauge- invariant regularization of the theory, 
one possibility being the use of a spacetime lattice; cf. Refs. [1, 2]. 

The perturbative quantum field theory based on the action (2.14) contains only tree and 
one-loop diagrams because the gauge field does not propagate and because there are no 
fermion self-interactions present. The effective gauge field action from the n = sector of 
the field theory (2.14) over the spacetime M 3 is directly related to the effective action from 
charged massless Dirac fermions over M 3 . 

The n = sector of the theory (2.14) has therefore the same "parity anomaly" as the 
standard M 3 theory [15-17]. The anomaly manifests itself in a contribution to the effective 
action T [B'} of the form 

fd 3 ys 7ruj C s[B' ,B[,B' 2 }, (2.15) 
Jm 3 

where u>cs is the Chern-Simons density 

u cs [B ,B 1 ,B 2 ] = J-, e K ^tr (B kX B, - lB K B x B^) , (2.16) 

in terms of the Yang-Mills field strength B K \ = d K B\ — d\B K + [B K , B\] , with indices running 
over 0, 1, 2, and the Levi-Civita symbol c kAm normalized by e 012 = +1 . The factor s in the 
anomalous term (2.15) is an odd integer which depends on the ultraviolet regularization 
used and we take s = +1 . Note that the Chern-Simons integral (2.15) is a topological 
term, i.e., a term which is independent of the metric on M3 . The total contribution of the 
n^O sectors to the effective action cannot be evaluated easily but is expected to lead to 
no further anomalies; cf. Refs. [1, 2]. 
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The anomalous term (2.15) gives a contribution to the four-dimensional effective action 
F[B'] of the form 

f d 4 x^ \-u C s[BiB[,B' 3 ] + -u;cs[B' ,B' 2 ,B' 3 ]\ , (2.17) 
J M 2np { p p J 

in terms of the usual Cartesian coordinates x^ and the corresponding four-dimensional gauge 
fields B' (x), p = 0, 1,2,3, and with the definition p 2 = (x 1 ) 2 + (x 2 ) 2 . In the form written, 
Eq. (2.17) has the same structure as the CPT anomaly term (4.1) of Ref. [1] for the R 3 x S 1 
manifold. The main properties of the anomalous term will be recalled in the next subsection. 

For completeness, we mention that the CPT anomaly also occurs for another type of 
orientable manifold, ^ = 1 x iV 3 , which is both multiply connected and geodesically com- 
plete. This particular space manifold A3 has a single wormhole, constructed from Euclidean 
three-space 1R 3 by removing the interior of two identical balls and properly identifying their 
surfaces (in this case, without time shift); cf. Refs. [11, 14]. The space "defect" is point- 
like if the removed balls are infinitesimally small. For an appropriate class of gauge fields 
B"(x), the parity anomaly gives directly an action term analogous to Eq. (2.17), now over 
the spacetime manifold iV and with a dipole structure for the uics terms in the integrand. 
The detailed form of this anomalous term is, however, somewhat involved (see Appendix 
A 1) and, for the rest of this section, we return to the original manifold M = M 4 \ M 2 . 

C. Abelian anomalous term 

The gauge field B^(x) will now be restricted to the Abelian Lie subalgebra u(l) C so(10) 
which corresponds to electromagnetism and the resulting real gauge field will be denoted 
by A^(x). In this case, the trilinear term of the Chern-Simons density (2.16) vanishes. 
For the ^-independent gauge fields (2.9) restricted to the u(l) subalgebra, the anomalous 
contribution (2.17) to the effective gauge field action T[A'] is simply 

__L f d 4 * e^ A g ; d v A' K {x) A' x (x) , (2.18) 
with the Levi-Civita symbol e fll ' KX normalized by e 0123 = +1, 

«>M-' = i (0,-J, J,0), (2.19) 

and p 2 = (a; 1 ) 2 + (x 2 ) 2 . 

Four remarks on the result (2.18) are in order. First, the action term (2.18), with q'^ = 
d^f, is invariant under a four-dimensional Abelian gauge transformation, A' (x) — > A'^ (x) + 
dn£(x). Second, the Lorentz and time-reversal (T) invariances are broken (as is the CPT 
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invariance), because the q'^ components in the effective action term (2.18) are fixed once and 
for all to the values (2.19); cf. Ref. [5]. Third, we expect no problems with unitarity and 
causality for the photon field, because the component q' vanishes exactly; cf. Refs. [18, 19]. 
Fourth, the overall sign of expression (2.18) can be changed by reversing the direction of the 
^-axis; cf. Eq. (2.15). 

After a partial integration, the anomalous contribution (2.18) can be generalized to the 
following term in the effective action T[A]: 



with the field strength F^x) = d^A u (x) — d u A^(x) and the integration domain extended 
to R 4 , which is possible for smooth gauge fields A^(x). (See Refs. [20-22] for a related 
discussion in the context of axion electrodynamics.) The factor fu{x\ A] in the anomalous 
term (2.20) is both a function of the spacetime coordinates x^ (on which the partial derivative 
8^ acts to give q^) and a gauge- invariant functional of the gauge field A^(x). This functional 
dependence of fu involves, most likely, the gauge field holonomies, defined as hc[A] = 
exp (i § c dx^ Ap{x)) for an oriented closed curve C; cf. Sec. 4 of Ref. [1]. 

Note that the functional Jm{x', A] in Eq. (2.20) is defined over M 4 but carries the memory 
of the original (multiply connected) manifold M, as indicated by the suffix. [The same struc- 
ture (2.20) has also been found for a manifold with a single static wormhole; see Appendix 
A 1.] Moreover, the absolute value of f M is of the order of the fine-structure constant, 



with the electromagnetic coupling constant e oc g . The general expression for fu{x\A] is 
not known, but /m can be calculated on a case by case basis [that is, the function f'(x) for 
the gauge field configuration A', f"(x) for A", et cetera]. 

III. MODELS 
A. Motivation 

The exact calculation of the CPT anomaly is prohibitively difficult for two or more 
punctures or wormholes. However, we expect the Abelian anomalies from the individual 
defects (considered in Sec. II and Appendix A 1) to add up incoherently, at least over large 
enough scales. We, therefore, assume that the total anomalous effect of the defects can be 
described by a contribution to the effective gauge field action T[A] of the form (2.20) but 
with f M (x;A] replaced by a background field gi(x). This background field gi(x) carries 




(2.20) 



\f M (x;A]\ = 0(a), a = e 2 /(An), 



(2.21) 
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the imprint of the topologically nontrivial structure of spacetime as probed by the chiral 
fermions. [The original spacetime manifold may, of course, have additional structure which 
does not contribute to the CPT anomaly and does not show up in gi(x).] 

As discussed in the Introduction, we consider in this paper a particular type of spacetime 
foam for which the defects of three-space are static and have randomly distributed positions 
and orientations. The average distance between the defects can be assumed to be small 
compared to the relevant scales (set by the photon wavelength, for example) and the detailed 
form of gi(x) is not important for macroscopic considerations. We, therefore, consider 
gi(x) to be a "random" field and only assume some simple "statistical" properties. These 
statistical properties will be specified in Sec. IV. 

It should, however, be clear that the background field gi(x) is not completely random. 
It contains, for example, the small-scale structure of the individual anomaly terms. The 
randomness of gi(x) traces back solely to the distribution of the static defects whose physical 
origin is unknown. In fact, the aim of the present paper is to establish and constrain some 
general characteristics of these hypothetical defects. 

In the rest of this section, we present two concrete models with random background fields. 
The first model describes the propagation of a single real scalar field and the second the 
behavior of the photon field. Both models are defined over Minkowski spacetime, M = R 4 
and g^(x) = rj^ . 

B. Real scalar field 

The scalar model is defined by the action 



where go is a real scalar background field of mass dimension zero. The background field go 
is assumed to be random and further properties will be discussed in Sec. IV. 
The corresponding equation of motion reads 



scalar 




(3.1) 



(□ + m 2 ) <j>(x) = -d^go(x) d^{x) , 



(3.2) 



with the following conventions: 



(O = diag(l, -1,-1,-1). 



(3.3) 
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There is no equation of motion for g , because go is a fixed background field [the random 
coupling constants go(x) are quenched variables]. The scalar field equation (3.2) will be seen 
to have the same basic structure as the one of the photon model in the next subsection. 
For the free real scalar field we define 

(O) (A;) = 27i5(k fI k >1 -m 2 ) [6(k°)ft°\k) 

+ 6(-k )ft \-k)*]. (3.4) 

The causal Green function A(x) of the Klein-Gordon operator, 

(n + m 2 )A(x) = 5 4 (x), (3.5) 

is given by 

a(x) = S Aik e ~*~ K[k) ■ 

with the Feynman prescription for the k° integration contour [23]. 
C. Photon field 

The photon model is defined by the action 

/photon =-\ \ d 4 X (F^(X)F^(X) 

+ 9l (x) F KX (x)F K \xj) , (3.7) 
where the Maxwell field strength tensor F^ v and its dual F kX are given by 

Fixv = d^A u — d u A fJi , 

F kX = \ e KX ^ , (3.8) 

with t KXhlu the Levi-Civita symbol normalized by e 0123 = +1 . 

The random (time-independent) background field g\ in the action (3.7) is supposed to 
mimic the anomalous effects of a multiply connected (static) spacetime foam, generalizing 
the result (2.20) for a single puncture or wormhole. Following Eq. (2.21), the amplitude of 
the random background field g± is assumed to be of order a. The typical length scale over 
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which gi varies will be denoted by Zf oam and further properties will be discussed in Sec. IV. 
Note that models of the type (3.7) have been considered before, but, to our knowledge, only 
for coupling constants varying smoothly over cosmological scales; cf. Refs. [18, 24]. 
The equations of motion corresponding to the action (3.7) are given by 

UA v {x) = -d^ 9l {x)F^{x), (3.9) 
provided the Lorentz gauge is used, 

& i A ll (x)=0. (3.10) 

The random coupling constants gi(x) in the action (3.7) are considered to be quenched 
variables and there is no equation of motion for g±. As mentioned above, the basic structure 
of Eq. (3.9) equals that of Eq. (3.2) for the scalar model. 
The general solution of the free field equation, UA^ = 0, is 

4 0) (*) = / d4 *^ e "^' 

Af(k) = 2n5(k,kn[e(k°)A^(k) 

+ e(-k°) If (-£)•], (3.ii) 

with 

Fl(°)(A;)=0, (3.12) 

due to the gauge condition (3.10). The Feynman propagator [23] for the free field is simply 
A^ v {x) = irj^y A(x), with A(x) given by Eq. (3.6) for m = . 

IV. RANDOM BACKGROUND FIELDS 

The two models of the previous section have random background fields go and g±. For 
simplicity, we assume the same basic properties for these two random background fields and 
denote g and g 1 collectively by g in this section. 

A. General properties 

The assumed properties of the background field g(x) are the following: 

1. g is time-independent, g = g(x), 

2. g is weak, \g(x) \ 1, 
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3. the average of g(x) vanishes in the large volume limit, 

4. g(x) varies over length scales which are small compared to the considered wavelengths 
of the scalar field and photon field , 

5. the autocorrelation function of g(x) is finite and isotropic, and drops off "fast enough" 
at large separations (see Sees. IV C and V below). 

The random background field gi(x) for the photon case is considered to incorporate the 
effects of a multiply connected, static spacetime foam (cf. Sec. Ill A) and assumption 5 about 
the lack of long-range correlations can perhaps be relaxed. Indeed, long-range correlations 
could arise from permanent or transient wormholes in spacetime; cf. Refs. [13, 14]. But, for 
simplicity, we keep the five assumptions as listed above. 



B. Example 

A specific class of random background fields g(x) can be generated by superimposing 
copies of a localized, square-integrable shape h(x) with random displacements. The dis- 
placements x n are uniformly distributed over a ball of radius R embedded in M 3 and have 
average separation a. The number of elementary shapes is then given by iV = (4ii/3)R 3 /a 3 . 

Concretely, we take 

N 

g N (x) = a ^ e n h(x - x n ) , (4.1) 

n=l 

where the numbers e n = ±1 are chosen randomly, so that 5^(5) vanishes on average. The 
background field g(x) is defined as the infinite volume limit, 

g(x) = lim g N (x) , (4.2) 

where a is held constant. 

The mean value and autocorrelation function of g are defined by 

(g) = lim - 1 / d 3 x g(x) , (4.3) 



R^oo (4tt/3)R 3 j 

\x\<R 



S A3gm9[s+l) - (4 - 4) 



]<R 



In the limit R — > oo, we have perfect statistics (N — > oo) and find 

(g) = 0, (4.5) 
C(x) = a 2 a' 3 [ d 3 y h(y) h(y + x) . (4.6) 

JR3 
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The Fourier transforms of h and g^ are given by 



h(k) = I d s xe- ik - £ h(x), 



(4.7a) 



N 



9n{ 



n=l 



(jfc) = a^e n / d 3 f e- ife "/i(f -f n ) 



= aVN h(k)G N {k), 



(4.7b) 



with 



JV 



v n=l 



e n e 



(4.8) 



The function (4.8) can be decomposed into an absolute value and a phase factor. The 
absolute value fluctuates around 1, as follows from the expression 



G N {k) 



6777, 6 



\k-(x n Xm) 



(4.9) 



where the double sum on the right-hand side scatters around 0. The fluctuation scale of 
|Gjv(/c)| is of order 2n/R and the same holds for the phase of Gat(/c); see Fig. 1. 



C. Random phase assumption 

For a given random background field g(x) over IR 3 , we define the truncated Fourier trans- 
form g R (k) by 

g R {k)= j d 3 xe- ils g(x), (4.10) 

|x|<fl 

with an implicit dependence on the origin of the sphere chosen. This function g R can be 
parametrized as follows: 

g R (k) = \/NP H(k) G R (k) , (4.11) 
where the real function H{k) is obtained from the finite autocorrelation function (4.4), 

(2^)3 jd 3 k^ s H(k) 2 , (4-12) 

with the sign of H{k) chosen to make the function as smooth as possible (i.e., without cusps). 
Because the background field g(x) is real, we have also H(—k) G R {—k) = H{k) G R (k)* . 

- — ■ — * 

Motivated by the results of the preceding subsection, we assume G R (k) to vary over 
momentum intervals of the order of 2%/R and \G R \ to have unit mean value. Furthermore, 
H(k) 2 is taken to be smooth on scales of the order of 2n/R, provided R is large enough. For 
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FIG. 1: The top and middle panels show the momentum-space function Gjv(fe), defined by the 
one-dimensional version of Eq. (4.8), for N = 100 and R = 500. The bottom panel shows the 
corresponding position-space function gw(x), defined by the one-dimensional version of Eq. (4.1), 
for the profile function h(x) = (20/9) (sin x/2 — sin x/20) jx. The Fourier transform h{k) of this 
particular profile function is nonzero and constant for \k\ G [fci ow , fchigh], with fchigh = 10 k\ ow > 0, 
and zero otherwise. 
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the isotropic case H{k) = H(k), with k = \k\, and large enough R, we explicitly assume the 
inequality 

\ dk cos(kR)H(k) 2 </ dkH(k) 2 , (4.13) 
Jo Jo 

with a finite positive number on the right-hand side. 

Now suppose that we have to evaluate a double integral of the form 



1^1 



d 3 pf(p,k)G R (k)G R (p), (4.14) 



where / is a function which is approximately constant over momentum intervals of the order 
of 2n/R. Due to the rapid phase oscillations of Gr, a significant contribution can arise only 
for k ~ — p and G R (k) G R (p) can effectively be replaced by a smeared delta-function. 

The specific form of this smeared delta-function does not matter in the limit R — > oo and 
we simply choose 



j=l 3 

We then have the result 

Jd 3 k Jd 3 pf(p,k)G R (k)G R (p) 



Jd 3 k jd 3 pf(p,k)5 3 R (k + p), (4.16) 

with a normalization factor J\f to be determined shortly. On the right-hand side of Eq. (4.16), 
the absolute values of G R have been replaced by their average value 1, since f(p, k) is assumed 
to be slowly varying. 

In order to determine the normalization factor A/", reconsider the autocorrelation function 
(4.4). Using Eqs. (4.11) and (4.15), the autocorrelation function can be written as follows: 

C(x) = lim — i- / d 3 kd 3 p 5 3 R (k + p) 
R—>oc (27r) ,i J 

xH(k)H(p) G R (k) Gr(p) ■ (4.17) 

Since there is already a delta-function present in the integrand of (4.17), the product 
G R (k)G R {p) can simply be replaced by M5 3 R (k + p), which gives the correct behavior for 

— * 

k ~ —p. With 

4(£ + V) S R (k + p)->?j5 3 (k + p), (4. 18) 

7r° 

for R — > oo, Eqs. (4.12) and (4.17) are consistent if 

A/" = (vr/i?) 3 . (4.19) 

This fixes the normalization factor Af of Eq. (4.16), which will be used extensively in the 
next section. 
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V. DISPERSION LAWS 
A. Scalar model 

We first turn to the scalar model (3.1), as a preliminary for the calculation of the photon 
model in the next subsection. In order to have a well defined Fourier transform of go, the 
system is put inside a sphere of radius R with free boundary conditions for the real scalar 
field 0. In momentum space, the scalar field equation (3.2) becomes 



where Gqr is assumed to have the statistics properties discussed in Sec. IV C. Furthermore, 
we have assumed isotropy of H ; see Sec. IV A. 

We will now show that the main effect of the random background field go(x) can be 
expressed in the form of a modified dispersion law for the scalar modes, at least for large 
enough wavelengths. The basic idea is to expand the solution of Eq. (5.1) perturbatively 
to second order in g (there is no contribution at first order) and to compare it with the 
first-order solution of a modified field equation, 




(5.1) 



According to Eq. (4.11), the momentum-space function g 0R has the form 



g 0R (k) = 27n5(A; ) H (\k\) G 0R (k) , 



(5.2) 



(k 2 - m 2 ) <f)(k) = A(\k\) (j)(k). 



(5.3) 



The dispersion law will then be given by 



u; 2 = m 2 + \k\ 2 + A(\k\) , 



(5.4) 



with A expressed in terms of the random background field go(x). 



In the limit R — ■> oo, the perturbative second-order contribution to 4>{k) reads 




P), (5.5) 
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where Eq. (5.2) has been used. According to Eqs. (4.16) and (4.19), we can replace the 
product GqrGqr by (n/R) 3 5 R and take the limit R — > oo, 

(2) W = I ^p A(A;) Jd 3 q H {\q\fMk\k-q) (q (q '- kj) (q ■ k) j<® '(*) . (5.6) 

At the given perturbative order, k^k^—m 2 can be put to zero in the integrand of Eq. (5.6), 
because of the delta-function contained in the free field <j>o(k) of Eq. (3.4). By comparison 
with the first order solution of Eq. (5.3), the operator A is identified as 

Following the discussion of Sec. IV A, we now assume that H (\q\) vanishes for momenta 

— * 

\q\ < Qiow (cf. Fig. 1) and that \k\ < q\ ow /2. Performing the angular integrals in Eq. (5.7), 
one finds 

q-2\k\ 



g + 2|A;| 



(5.8) 



where the lower limit of the integral is effectively q\ ow . Because of the large- wavelength 
assumption \k\ < q\ ow /2, the same result is obtained if the causal Green function A in the 
first integral of Eq. (5.5) is replaced by, for example, the retarded or the advanced Green 
function. 

— # 

Next, the logarithm of Eq. (5.8) is expanded in terms of \k\/q, which gives 

^(1*1) = 2^ Ho{qf [~\ q 2 \k\ 2 ~\ \k\^j + O (|^| 6 /ft 4 ow) • (5-9) 

This result can also be written as 

A ^ = ~ Ys Co(0) ^ ~ if ( [ dx x Co{x) ) + ° ^ /qL ) ' (5 - 10) 

in terms of the autocorrelation function (4.12), which is isotropic because H is, 

C = C (|f|) = j£y 2 fdq q* ^ Mqf . (5.11) 

In order to make the identification of the |/c| 4 prefactor in Eq. (5.10), it has been assumed 
that the behavior of H (q) 2 is sufficiently smooth; cf. Eq. (4.13). 
From Eqs. (5.4) and (5.10), the dispersion law of the scalar is 

to 2 = m 2 + (1 - oo) \k\ 2 - b \k\ 4 + ... , (5.12) 

where the dots stand for higher order contributions and the positive coefficients a and 6 
are given by 

7r 2ir f°° 

a o = -C o (0), b o = j^J o dxxC (x). (5.13) 
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The isotropic autocorrelation function C (x) is defined by Eq. (4.4) with g replaced by the 
random background field go from the action (3.1). 

The main purpose of the scalar calculation is to prepare the way for the photon calculation 
in the next subsection, but let us comment briefly on the result found. Setting the scalar 
mass m to zero in the action (3.1), our calculation gives no additional mass term in the 
dispersion law (5.12). It is, however, not clear what the random background field g of 
the model (3.1) really has to do with a (static) spacetime foam. The propagation of an 
initially massless scalar could very well be strongly modified in a genuine spacetime foam; 
cf. Ref. [10]. 

B. Photon model 

To obtain the dispersion law for the photon model (3.7), we proceed along the same lines 
as for the scalar case. Again, the system is put inside a sphere of radius R, so that the 
truncated Fourier transform g iR occurs in the momentum-space field equation, 

k 2 A»(k) = - JL J d 4 g g m {q) e^q, (k - q) K A x (k - q) . (5.14) 
For this truncated background field g 1R , we assume a form analogous to Eq. (5.2), 

g 1R (k) = 2nS(k°) y/(An/3)R 3 H^k]) G 1R (k) , (5.15) 

where G± R is a random function of the type discussed in Sec. IV C. 

The second-order contribution to the perturbative solution of Eq. (5.14) is given by 

AP»{k)= lim 7^A(A;) /d 4 g g 1R (q) q, (k - q) K A(k - q) 

X Jd 4 P g 1R (k- P -q)e a(37X (k- P -q) a P (S A^( P ). (5.16) 

Using Eq. (5.15) and replacing the product Gi R G\ R by (n/R) 3 S R , one finds in the limit 
R — > oo: 

A {2)u (k) = -A(k) B\(k) A^\k) , (5.17) 

with 

B ^*)= 12^Fp I^HM? (fc _ g )2 + ie W] <b(k-q) K <fk0, (5.18) 

where q° = and the square brackets around the indices af3^ denote antisymmetrization 
with unit weight. 



18 



Since A^ 1 (k) in Eq. (5.17) contains a factor S(k lx k^ i ) and furthermore obeys the Lorentz 
gauge condition (3.12), we have 



B\ A {0 > = , 5% = B n m A {0)m , 



with 

B n m (k) 

and 
B{\k\) 



12 (2nf 



dV^i(kl) 2 



2A;-g - |g] 2 + ie 



l&l 2 



(5.19) 
(5.20) 



1 



24tt Jo 



dq q 2 H x {qf l-f 



8\k\ 



In 



q-2\k\ 



q + 2\k\ 



(5.21) 



where the indices m and n run over 1, 2, 3. In order to perform the angular integrals in 
Eq. (5.20), we have again assumed that Hi(\q\) vanishes for momenta \q\ < q\ ow (cf. Fig. 1) 
and that \k\ < q\ ow /2. Making an expansion in \k\/q for the logarithm in Eq. (5.21), we 
obtain 



B(\k\) 



dq Hi(qf 



2 I7J2 



q z \k 



24vr J 

which is equivalent to the scalar result (5.9). 

The modified field equations in momentum space are thus 



\k\ A + ... 



(5.22) 



k 2 A°(k) = 0, k 2 A m (k) = B m n (k) A n {k) . (5.23) 

It follows that the dispersion law for the scalar and longitudinal modes remains unchanged, 

uj 2 = \k\ 2 , (5.24) 

whereas the one for the transverse modes is modified, 

u 2 = \k\ 2 + B(\k\). (5.25) 

Of course, only the transverse modes are physical; see, e.g., Sec. 3.2 of Ref. [23]. From Eqs. 
(5.22) and (5.25), the dispersion law for the (transverse) photons is 



(1 - a,) c 2 \k{ 2 -b lC 2 \k\ A + 



(5.26) 



where the bare light velocity c has been restored and the positive coefficients ai and bi are 
given by 

7T 2-7T f°° 

ai = — d(0), h = — dxxC^x). (5.27) 
18 ID J 
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The isotropic autocorrelation function C±(x) is defined by Eq. (4.4) with g replaced by the 
random coupling constant g\ from the action (3.7). 

According to the discussion in Sees. II C and III C, the random background field g\{x) 
has a typical amplitude of the order of the fine-structure constant a and a typical length 
scale of the order of lf oam . The coefficients (5.27) may therefore be written as 

ax = a 2 7x , &i = a 2 l 2 oam , (5.28) 

with a positive constant 71. This last equation defines, in fact, the length scale lf 03ja . For 
the example background field of Sec. IV B with isotropic profile function h = h(\x\), one 
has /foam oc lh(lh/a) 3 ^ 2 , where lh is the typical length scale of the profile function h; cf. 
Eq. (4.6). For this type of background field, the ratio lf oam /a can be large if lh ^> a . In 
Appendix A 2, another example background field is discussed, which is especially tailored 
to the case of permanent wormholes (average distance a between the different wormholes, 
effective transverse width 2 lh for the individual wormholes, and long distance d between 
the individual wormhole mouths). The result for Zf oam is again found to be proportional to 
lh (lh/a) 3 ^ 2 and only weakly dependent on the individual scale d. 

The dispersion law (5.26) violates Lorentz invariance but not CPT invariance. There is 
no birefringence, as the background field gi(x) is assumed to have no preferred direction. 
Note also that motion of the detector relative to the preferred frame defined implicitly by 
the static background gi(x) would bring in some anisotropy but, at the order considered, 
no birefringence (the two modes propagate identically in the instantaneous rest frame of the 
detector). 

The modifications of the photon dispersion law found in Eq. (5.26) are rather mild (see 

— * 

Ref. [25] for some general remarks on the absence of odd powers of | A: J in modified dispersion 
laws). These modifications are not unlike those of Ref. [10] for a simply connected, nonstatic 
spacetime made-up of many S 2 x S 2 "gravitational bubbles." More drastic changes in 
the photon dispersion law have, for example, been found in certain loop quantum gravity 
calculations [26, 27]. Compared to these calculations (which have the Planck length as the 
fundamental scale), ours is relatively straightforward, the only prerequisite being a multiply 
connected topology which is then probed by the chiral fermions of the Standard Model; see 
Sec. VII for further discussion. 

VI. EXPERIMENTAL LIMIT 

In the previous section, we have derived the photon dispersion law for the simple model 
(3.7). Using the definitions (5.28) and considering a to be parametrically small, we have the 
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following expression for the quadratic and quartic terms in the dispersion law (5.26): 

^4n^-4n« 2 Cn^, (6-1) 

with the renormalized light velocity 

c rcn = c a/1 - a 2 7i (6.2) 

and the simplified notation k = \k\ . 

The group velocity v g (k) = du;/d/c is readily calculated from Eq. (6.1). The relative 
change of v g {k) between wave numbers k\ and k 2 is then found to be 



Ac 

c 



ki,k 2 



Vg(h) ~ V g (k 2 ) 



2 1^-^1 « 2 IL, (6.3) 



where Ac/c is a convenient short-hand notation and a pa 1/137 the fine-structure constant. 

As realized by Amelino-Camelia et al. [28], one can use the lack of time dispersion in 
gamma-ray bursts to get an upper bound on Ac/c. But for our purpose, it may be better 
to use a particular TeV gamma-ray flare of the active galaxy Markarian 421, as discussed 
by Biller et al. [29]. Schaefer [30] obtains from this event 

Mkn 421 

< 2.5 x 1(T 14 . (6.4) 



Ac 

c 



fci=2.5xl0 16 cm- 1 
fe 2 =1.0xl0 17 cm- 1 



The reader is invited to look at Fig. 2 of Ref. [29], which provides the key input for the 
bound (6.4), together with the galaxy distance D. In fact, the right-hand side of Eq. (6.4) 
is simply the ratio of the binning interval for the gamma-ray events (At pa 280 s ) over the 
inferred travel time (D/c pa 1.1 x 10 16 s). 

Combining our theoretical expression (6.3) and the astrophysical bound (6.4), we have 
the following "experimental" limit: 

/ foam < 1.6 x 10- 22 cm, (6.5) 

where / foam is defined by (5.27) and (5.28), in terms of the autocorrelation function (4.4) for 
the static random variable gi(x) from the action (3.7). In the next section, we comment 
briefly on the possible interpretation of this result. 

Since the spacetime foam considered in this paper has no preferred spatial direction, 
experimental limits obtained from bounds on the birefringence of electromagnetic waves 
[18, 31, 32] do not apply. [See also the remarks in the penultimate paragraph of Sec. VB.] 
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Note that the background field gi(x) vanishes on average and that the modification of the 
dispersion law in Eq. (5.26) is a second-order effect, governed by the autocorrelation of the 
background field. As shown by Eq. (6.1), the fundamental length scale /f oam appears only 
at quartic order in the photon wave number, which makes /f oam difficult to constrain (or 
determine) experimentally. 

VII. CONCLUSION 

The present article contains two main theoretical results. First, it has been shown in 
Sec. II and Appendix A 1 that two particular types of "defects" of a noncompact spacetime 
manifold generate an anomalous breaking of Lorentz and CPT invariance. The anomalous 
term of the effective gauge field action can be written in a simple form, Eqs. (2.20) and 
(2.21). The spacetime defect is thus found to affect the photon field far away, the agent 
being the second-quantized vacuum of the chiral fermions. Inversely, the CPT anomaly can 
be used as a probe of certain spacetime structures at the very smallest scales. 

Second, a modified photon dispersion law, Eq. (5.26), has been found in a model which 
generalizes the single-defect result. The random background field gi(x) of this model (3.7) 
traces back to a postulated time-independent foam-like structure of spacetime consisting of 
many randomly oriented and randomly distributed defects. The static random background 
field gi(x) of the photon model breaks Lorentz and CPT invariance and selects a class 
of preferred inertial frames. Incidentally, the calculated photon dispersion law shows the 
Lorentz noninvariance present at the microscopic scale but not the CPT violation. 

This article also gives an "experimental" result. Following up on the suggestion of Refs. 
[28, 29], it is possible to use observations of gamma-ray bursts and TeV flares in active 
galactic nuclei to obtain an upper bound (6.5) on the length scale /f oam of the random 
background field gi(x) of the model considered. As such, this upper bound constitutes a 
nontrivial result for a particular characteristic of a multiply connected space manifold at 
the very smallest scale, granting the relevance of the model (3.7) for the effects of the CPT 
anomaly (see, in particular, the discussion in Sec. Ill A). 

The upper bound (6.5) on /f oam is, of course, eleven orders of magnitude above the Planck 
length, /pi an ck = a/ Gh/c 3 pa 1.6 x 10 -33 cm. But it should be realized that we have no real 
understanding of the possible topologies of spacetime, be it at the very smallest scale or the 
very largest. It is even possible that Zf oam and Zpi anc k are unrelated; for example, if /f oam is 
not a quantum effect, but has some other, unknown, origin. 
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APPENDIX A: CALCULATIONS FOR STATIC WORMHOLES 
1. CPT anomaly from a single wormhole 

The orientable three-space N$ considered in this subsection has a single permanent worm- 
hole (traversable or not). The wormhole is constructed by removing two identical open balls 
from IR 3 and properly identifying their surfaces; see, e.g., Sec. II A of Ref. [11] and Sec. 
15.1 of Ref. [14] for further details. The removed balls have diameter b and are separated 
by a distance d in M 3 . The length of the wormhole "throat" is then zero, whereas the long 
distance between the wormhole "mouths" is d. As a particularly simple case, the width b of 
the wormhole mouths is put to zero, so that the three-space N$ is essentially R 3 with two 
points identified. The coordinates of these identified points are taken to be x = ±(d/2) f, 
for an arbitrary unit vector f. 

Like the punctured three-manifold M 3 of Sec. II, this three-space N 3 is multiply connected 
and the CPT anomaly [1] may be expected to occur. In order to show this, we proceed 
along the same lines as in Sec. II. That is, we introduce a suitable coordinate system over 
the spacetime manifold N = R x N 3 and restrict our attention to a class of background 
gauge fields which allows us to trace the CPT anomaly back to the three-dimensional parity 
anomaly [15-17]. 

Let $(x) be a scalar function on M 3 , so that the integral curves of V$ are noncontractible 
loops on N 3 . As a concrete example, take 



which resembles the potential of an electric dipole. The manifold N$ is now parametrized 
by the coordinates 77, p, 0, where G [0, 27r) is the azimuthal angle around the r-axis, 
77 G [— 7r/2,7r/2) is defined by 77 = $(x), and p G [0, 00) is a coordinate perpendicular to 
and 77. The pair (p, 0) parametrizes the "equipotential" surfaces of $. The equipotential 
surface has the topology of a two-sphere for 77 7^ and that of a plane for 77 = (p and 
now correspond to the usual polar coordinates); see Fig. 2. 

For our purpose, it suffices to establish the CPT anomaly for one particular class of 
gauge fields. We, therefore, take the gauge fields to be independent of 77 and to have no 
component in the direction of 77. These gauge fields will be indicated by a double prime. 
With appropriate vierbeins, the fermionic fields are assumed to have periodic boundary 
conditions in 77. 

The anomaly can be calculated for the plane 77 = by making use of the three-dimensional 




(Al) 
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FIG. 2: Parametrization of the three-space N$ with a permanent wormhole at x = ±1 and y = z = 
(x and y coordinates shown, z coordinate suppressed). The thin lines indicate "equipotential" 
surfaces on which p runs from to oo and (j) from to 2ir (shown here are <j) = and (j> = ir). The 
thick lines show noncontractible loops on which r\ runs from —ir/2 to +ir/2. The minimal length 
of these noncontractible loops is given by d in Eq. (Al), which has the value 2 here. 

parity anomaly. Choosing f = (0, 0, 1), the relevant Chern-Simons term reads 

^cs[-EC B", B'l] 

= ^ <l'l ^ tr {BIB: - | B'*B'Ib:) , (A2) 

with 

tf")= (o, e v | v=Q ) =(0,0,0,-1), (A3) 

where e,, denotes the unit vector in the 77-direction, e v = V$/|V<1?|. The three-dimensional 
anomaly term is then 

/ dx° dx 1 dx 2 Tr s ^cs [Bq, B'l, B'±] , (A4) 

with an odd integer s from the ultraviolet regularization (in the following, we take s = 1). 
With the vector 

<j \B"\ = ^ ^ tr (B; U B: - I B'IB" V B" K ) , (AS) 
the anomaly term (A4) may be written as 

dx° [ Q[B"\ ■ dS , (A6) 

where denotes the surface of constant rj. 
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Since the fields considered are independent of rj, the result (A6) is independent of the 
particular equipotential surface chosen. Hence, 

- fdx° ! lo[B"} ■ dS (A7) 

is independent of 77 and can be averaged over 77. The anomaly term can then be written as 
a four-dimensional integral over the spacetime manifold i\f = R x iV 3 , 



j dx° j Q[B"} ■ dS 

_I r /2 dr] fdx° [ u)[B"] ■ dS 
71 J — 7r/2 J Jy„ 



1 



tt/2 

d 4 x |V$(£)| {u[B"\ ■ e v ) 



N 



= ^ J N d 4 x $(f ) V • u[B"\ , (A8) 

where R[&] is the range of values of the function i.e., the length of the coordinate 

interval of rj. For the particular function $ given in Eq. (Al), one has _R[$] = ir. 

If B"^ contains only an electromagnetic component A" , the anomalous contribution to the 
effective action is 

S/ X ^ e ^ F "^W' ( A9 ) 
up to an overall constant of order 1, which depends on the details of the theory considered. 
Here, the field strength F''{x) is defined by d^A'^x) — d u A'^(x) and $(x) is the scalar 
function (Al), possibly with a different normalization. The general structure of the term 
(A9) is of the form as given by Eq. (2.20), now with a functional /W(a;; A] in the integrand. 



2. Random background field and if oa m 



In this subsection, we present a simple background field g\{x) to mimic the anomalous 
effects from a random distribution of wormholes. The photon model (3.7) incorporates, 
in this way, the basic features of the CPT anomaly for a single wormhole as found in the 
previous subsection of this appendix. For purely technical reasons, we use, instead of the 
function (Al), a function $ which is the direct difference of two "monopole" contributions, 



$(£) = K 



h[x — -r 



— h [x + - r 



(AlO) 



with a displacement parameter 8 > and a normalization factor K > (see below). In ad- 
dition, the profile function h(x) is assumed to be isotropic, monotonic, finite, and integrable. 
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The random background field is a superposition of these "dipole" potentials, with random 
locations (x n ) and directions (f n ). The random background field also includes a factor a; 
cf. Eq. (A9). 

For iV = (4/3)"7r.R 3 /a 3 dipoles in a ball of radius R (with an average separation a between 
the different wormholes), we have 



9n 



\x) = aK ^2 \h(x- x n - -r n ) 



n=l 



-h(x-x n + |f„) 

where x n are random vectors with \x n \ < R and f n randomly chosen unit vectors. The 
final background field for the photon model (3.7) is given by the infinite volume limit, 
gi(x) = liniAr^oo Qn(x). The relevant autocorrelation function is 

1 



(AH) 



Ci( 



(4/3)7^3 J lsl<R 
,K 2 



d x g N (x) g N (x + y) 



d 3 x h^x — -fj — h(^x + -fj 



x 



h(x + y — - r j — h(^x + y + ^f^j 



(A12) 



where av^ denotes the average over the orientations r. 

The typical length scale which enters the photon dispersion law (5.26) is defined by Eqs. 
(5.27) and (5.28), 

1 



;2 

foam ~ 30a 2 



After a suitable shift of integration variables, one finds 

1 (A-B), 



I 2 

foam 
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in terms of two integrals, 

A = K 2 [ d 3 yd 3 x —r\h(x)h(x + y) , 

J \y\ a 

B = K 2 J d 3 yd 3 x 



1 



h(x)h(x + y) . 



(A13) 

(A14) 

(A15a) 
(A15b) 



y — f 5\ a 3 

The integral B is, in fact, independent of the direction f and there is no need to average over 
it. By using spherical coordinates (y, 9 y , <p y ) with respect to the r-axis for y and (x, 9 X , (f> x ) 
around the y-axis for x, one finds after performing three angular integrals: 



B 



A + K 



2 2 (27T) 2 



2„,2 



dy dx d9 x sin 9 X x y 



x ( ^ - i ) h(x) h(x + y) 



(A16) 
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FIG. 3: Long distance d between the wormhole mouths vs. parameter 5 entering the profile functions 
h of Eqs. (A10) and (A18), both in units of the width l^. 

With Eq. (A 14), this gives 



In order to get a concrete result for Zf oam , we choose the following profile function: 



with a scale parameter lh > 0. It is important to note that the peaks of $(x) from Eq. (A10) 
are not located at ±(5/2) f, but rather at ±(d/2) f for an effective distance d = d(S, lh)', see 
Fig. 3. The normalization factor K = K(5,lh) in Eq. (A10) is chosen such that <$>(x) has 
values ±1 at these peaks. The expressions for d and K are somewhat involved and need not 
be given explicitly. 

For the profile function (A18), we finally obtain 



with average separation a between the different wormholes and parameters 5 and lh for the 
individual wormholes; cf. Eqs. (A10) and (A18). As mentioned above, precisely the quantity 
Zf oam enters the quartic term of the photon dispersion law (6.1). For 5 /lh — > oo, one finds 





(A17) 




(A18) 




(A19) 



foam 




(A20) 
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FIG. 4: Length scale l{ oam entering the photon dispersion law (6.1) vs. individual wormhole scale 
d, with the further definition fi 2 = (tt 2 /15) l\/o?. 

which is essentially the same result as for a random distribution of "monopoles" with profile 
function h(x). For 5/lh j 0, on the other hand, the effective distance d approaches the value 
V / 4/5 l h and lf oam becomes 



44^=^1 x (0.279)*. (A2!) 



Figure 4 shows the rather weak dependence of /f oam on the individual wormhole scale d. 
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